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Abstract. The experimental data on F2 structure functions of the proton and deuteron were used 
to construct their moments. In particular, recent measurements performed with CLAS detector at 
Jefferson Lab allowed to extend our knowledge of structure functions in the large-x region. The phe- 
nomenological analysis of these experimental moments in terms of the Operator Product Expansion 
permitted to separate the leading and higher twist contributions. Applying nuclear corrections to 
extracted deuteron moments we obtained the contribution of the neutron. Combining leading twist 
moments of the neutron and proton we found d/u ratio at jc — > 1 approaching 0, although 1/5 value 
could not be excluded. The twist expansion analysis suggests that the contamination of higher twists 
influences the extraction of the d/u ratio at x — > 1 even at Q 2 -scale as large as 12 (GeV/c) 2 . 
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DATA ANALYSIS 

QCD through the Operator Product Expansion (OPE) allows to relate measurable mo- 
ments of nucleon structure functions to the series of local operators, so called twists. 
The Leading Twist (LT) (first term in the series) represents the asymptotic freedom do- 
main. This term is completely determined by perturbative calculations and Lattice sim- 
ulations 1 . Higher Twists (HT) (all further terms) describe the virtual photon scattering 
off interacting partons. The complexity of this interaction and therefore of correspond- 
ing QCD operators increases with twist order. Calculations of these terms have been 
performed only in a few cases. 

The experimental data on the structure function moments were obtained recently for 
the proton and deuteron in Refs. [d and |0], respectively. Moments were extracted 
from a combined analysis of the structure functions F% measured at CLAS and other 
world data on the inclusive lepton-nucleon scattering. This was performed by integrating 
experimental data points independently at each fixed Q 2 value. Hence, the obtained Q 2 - 
evolution of the moments is not affected by any model dependence and can be directly 
compared to pQCD predictions. Example of measured moments is shown in Fig. [2 As 



Lattice simulations up to now are limited to a few lower moments 
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FIGURE 1. Total experimental n = 6 moments of the proton (full triangles) and deuteron per nucleon 
(open squares) structure functions F%, 



one can see the proton and deuteron moments have similar <2 2 -behavior, but different 
absolute value. Indeed, in the QCD at the LT two moments should have the same Q 2 - 
evolution. This observation suggests that the duality is valid for both targets and the 
contribution of nuclear HTs 2 is small in the covered <2 2 -range. 

We analyzed these experimentally extracted moments of the proton and deuteron 
structure functions F2 OLIjj to separate LT and HT terms. This was performed by fitting 
the data with the following expression: 

M n (Q 2 )=J Q l dx^- 2 F 2 (x,Q 2 )=LUa s ) + trt^ (^Y , (1) 

where as is the running coupling constant, /1 2 is an arbitrary scale (taken to be 10 
(GeV/c) 2 ), a x n is the matrix element of corresponding QCD operators, y% is the anoma- 
lous dimension, x is the order of the twist and k is the maximum HT order considered. 
The number of HT terms (k) in the expansion is of course arbitrary because we don't 
know at which l/<2 2 power the series converges. This prevents an evaluation of each 
separate HT term from the data. Instead, the total contribution of HTs can be extracted 
with good precision. In Fig. El one can see that taking two, three or four HT terms in 
Eq.[l]does not change the total HT contribution. Moreover, this result was expected be- 
cause the total HT contribution represents simply the difference between the data and 
calculated LT. 

Extracted LT components of the proton and deuteron moments can be combined 
now to obtain moments of the neutron structure function F%. In the Euclidean space 
of moments the convolution of the nuclear Impulse Approximation (IA) transforms into 



nuclear higher twists are mostly related to Final State Interactions of the nucleon in the nucleus 



-2 



xlO 



0.12 



£ 0.06 - 



0.02 - 



0.04 - 



0.08 - 



0.1 - 








10 



10 



Q 2 (GeV/c) 2 



FIGURE 2. Total higher twist contribution in n = 8 moment obtained with the present procedure for 
different number of terms in the OPE series (see Eq.[Q: solid line - two HT terms, dashed line - three HT 
terms, dotted line - four HT terms. 

a product of moments. This allows for a simple extraction of neutron moments from the 
following algebraic relation: 



where M„, M% and M„ are moments of the proton, neutron and deuteron, respectively. 
N„ is the moment of the nuclear momentum distribution f D i.e. the structure function 
of the deuteron composed of point-like nucleons (see Ref. Oy for details). This nuclear 
structure function f D was obtained from the data on the deuteron wave function. 

The data on proton and neutron moments can be used to study the contribution of 
u and d quarks in the proton at x — > 1. Assuming that u and d quark distributions in 
the proton and neutron are the same 3 the d/u ratio at the leading twist accuracy can be 
related to the ratio of neutron to proton structure functions F^/F^ix — ► 1). This ratio, 
in turn, is equal to the ratio of moments of these structure functions M%/M%(n — > °°). 
The latter equality requires only that F/' n (;c — > 1) — > a p ' n (l —x) , which follows from 
the analyticity of the forward Compton amplitude in OPE. Instead, the existence of a 
finite F^/F^ix — ► 1) limit guarantees that exponents b p and b n are exactly the same. We 
constructed these ratios of moments and results are shown in Fig.|3j As one can see our 
data tend to the "standard" 1/4 value used in most of parton distribution fits. This value 
corresponds to the vanishing d/u ratio. However, the precision of our data does not allow 
to exclude 3/7 value corresponding to d/u ratio 1/5. In the Fig. |3] one also can see the 
impact of the HT contribution on the ratio at Q 2 = 12 (GeV/c) 2 . The previous analysis |0] 
performed in the x-space showed that applying similar nuclear corrections the ratio of 



though the number of u and d quarks is of course different 
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FIGURE 3. Ratio of the neutron to proton moments as a function of n at Q 2 = 12 (GeV/c) 2 : full triangles 
show the ratio at the leading twist obtained in Ref. [3], open circles represent the ratio of measured inelastic 
moments including the higher twist contribution. 

structure functions F^/F^x—^ 1) goes to 3/7 at Q 2 = 12 (GeV/c) 2 . No corrections on 
the possible HT contamination has been applied in this article, assuming that at such 
large Q 2 they are negligible. Therefore, if instead of LT part we take measured inelastic 
moments, containing also HT terms, and construct the same ratios we should confirm the 
result from Ref. |4|]. Indeed, in the Fig.|3]one can see that the ratio of moments including 
HTs tends to 3/7 at largest n (n = 12 corresponds to x values about 0.75). This result 
could also be deduced from the isospin independence of HTs observed in the Ref. |0]. 

CONCLUSIONS 

We obtained the experimental data on moments of proton and deuteron structure func- 
tion F%. In these data contributions of the leading and higher twists were separated via 
OPE analysis. By combining the proton and deuteron moments and applying nuclear 
corrections we extracted moments of the neutron structure function F<i. The ratio of the 
neutron and proton moments at large n is related to the ratio of d and u quark contribu- 
tions in the proton at large-x The obtained ratio is consistent with the asymptotic limit 
d/u^0dXx^ 1, which originates from the dominance of soft, non-perturbative physics 
at large-x. Nevertheless the alternative value of 1/5, derived from helicity conservation 
arguments, is not excluded. The HT contamination to the ratio, if not subtracted as in the 
present analysis, influences significantly the extracted d/u ratio also at large <2 2 -values. 

REFERENCES 

1. M. Osipenko, et al., Phys. Rev. D67, 092001 (2003). 

2. M. Osipenko, et al., Phys. Rev. C73, 045205 (2006). 

3. M. Osipenko, et al., Nucl. Phys. A766, 142 (2006). 



i 



0- 



3= 



3 



3/7 



1/4 



4. W. Melnitchouk, and A. Thomas, Acta Phys. Polon. R27, 1407 (1996). 



